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Abstract. In this paper, we investigate the long-time structure of the heat kernel on 
a Riemannian manifold M which is asymptotically conic near infinity. Using geometric 
microlocal analysis and building on results of Guillarmou and Hassell in |GHlj . we give 
a complete description of the asymptotic structure of the heat kernel in all spatial and 
^vq . temporal regimes. We apply this structure to define and investigate a renormalized zeta 

function and determinant of the Laplacian on M. 
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OV 

1. Introduction 

We study the heat kernel on asymptotically conic manifolds. Asymptotically conic mani- 
folds should be thought of as those complete manifolds which are approximately conic near 
^ ■ infinity. More specifically, we have the following definition [GHlJ: 

Definition. Let [M, g) be a complete Riemannian manifold without boundary of dimension 
n, and let M be the usual radial compactification of M. Let (N, h Q ) be a closed Riemannian 
manifold of dimension n — 1. We say that (M,g) is asymptotically conic with cross-section 
1 (N,h ) if in a neighborhood of dM , M is isometric to [0,S) X x N y with the metric 

oo : 

O- dx 2 h(x) 

9 = ^T + —- (!) 

Here x is a smooth function on M with x = and dx on dM (we call this a boundary 

(N 



defining function for dM) and a smooth family of metrics h(x) on N with h(0) = ho. 



In particular, Euclidean space MJ 1 is asymptotically conic with cross-section S n ; we may 
choose x = r -1 . Any complete manifold which is exactly Euclidean or conic near infinity is, 
of course, also asymptotically conic. The condition ([T|) may be weakened by replacing h{x) 
with any symmetric 2-tensor h'(x,y) which restricts to a metric ho{x) on the boundary at 
x = 0; an observation of Melrose and Wunsch [MW] shows that these conditions are in fact 
equivalent. 

Asymptotically conic manifolds are a relatively well-behaved class of manifolds, and as 
such the theory of the heat equation is relatively advanced. In particular, it is easy to see 
from ([T]) that all sectional curvatures of (M, g) approach zero as x goes to zero, and thus 
that (M, g) has bounded sectional curvature. For complete manifolds of bounded sectional 
curvature, a classical theorem of Cheng-Li- Yau |CLY] gives the following Gaussian upper 
bound for the heat kernel: 

Theorem 1. |GLYj There are nonzero constants C% and C<i such that the heat kernel on M , 
denoted H M (t, z, z'), satisfies 

H M (t,z,z')<^- 2 e-^. (2) 
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However, for many applications to spectral theory, one needs finer information about the 
structure of the heat kernel. The example we have in mind is the definition of the zeta 
function. Recall that if M is compact, the zeta function is defined for > n/2 by: 

1 f°° 

Cm(s) = I>)/ - I)** -1 df- (3) 

The zeta function has a well-known meromorphic continuation to all of C with a regular 
value at s = 0; the key is that the trace of the heat kernel has a short-time asymptotic 
expansion, which, along with the long-time exponential decay, enables us to write down an 
explicit meromorphic continuation [R]. The determinant of the Laplacian is then given by 
exp(— Ca/(0)); the determinant plays a key role in many problems in spectral theory, including 
the isospectral compactness results of Osgood, Phillips, and Sarnak |OPSR IOPS2t IOPS3] . 

We would like to define such a zeta function and determinant when M is asymptotically 
conic, with an eye towards applying these concepts to the spectral and scattering theory 
of asymptotically conic manifolds. There are several obstacles. First, the heat kernel is no 
longer trace class, so TrH M (t) does not make sense. Instead, we define the renormalized 
heat trace R TrH M (t) to be the finite part at 5 = of the divergent asymptotic expansion in 
5 of 

/ H M (t,x,y,x,y)x- n - 1 dxdy. (4) 

Jx>8 

Details, including the existence of this divergent asymptotic expansion, may be found in 
Section 3. We then formally define the renormalized zeta function: 

1 f°° 

r (m(s) = — / TrH (t)t s ~ l dt. (5) 
r (s) Jo 

However, to make sense of this definition and obtain a meromorphic continuation, we still 
need to understand the behavior of the renormalized trace - and hence of the heat kernel 
itself - as t — > and t — > oo. In particular, we need asymptotics in both the short and long 
time regimes. 

The short-time behavior of the heat kernel on an asymptotically conic manifold is rel- 
atively well-understood. Short-time heat kernels may be analyzed using techniques from 
semiclassical analysis. In this approach, the goal is to develop a 'semiclassical functional 
calculus' containing the heat kernel, modeled on standard semiclassical techniques as de- 
veloped, for example, in |DS] . The key functional calculus for this purpose, at least in the 
asymptotically Euclidean setting, is the Weyl calculus of Hormander |Hoj . 

An alternate approach, and one more suited to analysis of the renormalized trace and 
determinant, is to use geometric micro local analysis to first construct the heat kernel and 
then analyze its fine structure. The techniques of geometric microlocal analysis were first 
developed by Melrose and Mendoza to study elliptic PDE on manifolds with asymptotically 
cylindrical ends |MeMej . They have been extended by many other mathematicians and play 
a key role in the modern analysis of linear PDE on singular and non-compact spaces. In 
particular, in |Me3] . Melrose discusses some aspects of spectral and scattering theory on 
asymptotically conic manifolds. Albin, in [ATE] , uses these methods to investigate the short- 
time heat kernel on a variety of complete spaces, including asymptotically conic manifolds. 
His work can be used to obtain the fine structure that we need for the short-time heat kernel. 

The long-time problem is trickier: in the asymptotically conic setting, we no longer have 
exponential decay of the heat kernel as t — > oo. Indeed, from the structure of the Euclidean 
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heat kernel and ([2]), we expect that the leading-order behavior of H M (t, z, z') as t — > oo will 
be Ct~ n / 2 , and the leading-order behavior of the renormalized heat trace may be even worse. 
This lack of decay means that the zeta function may not be well-defined a priori for any s. 
We may split ([S]) into two integrals by breaking it up at t = 1, but there is no obvious reason 
for the integral from t — 1 to oo to have a meromorphic continuation to all of C. In order 
to obtain such a meromorphic continuation, we need an asymptotic expansion for the heat 
kernel as t — > oo. Moreover, we must understand how this expansion interacts with the heat 
trace renormalization. 

1.1. Main results. We solve this problem by using the methods of geometric microlocal 
analysis to obtain a complete description of the asymptotic structure of the heat kernel on 
M in all spatial and temporal regimes. Our key result is the following: 

Theorem 2. Let M be asymptotically conic. For any n > 2, and for any fixed time T > 0, 
the heat kernel on M is polyhomogeneous conormal on M 2 sc for t > T, where w = t^ 1 ^ 2 . 
The leading orders at the boundary hypersurfaces are at least at sc and n at each of bf , 
rbo, lbo, and zf, with infinite- order decay at lb, rb, and bf. 

This theorem gives a complete description of the asymptotic structure of the heat ker- 
nel for long time; previously, only estimates such as Theorem |CLYj were known. The 
space M 2 sc is a manifold with corners which is a blow-up of (w, x, y, x', y')-space; we give 
its explicit definition in Section 2, and also briefly discuss blow-ups and polyhomogeneous 
conormal functions. These concepts were originally introduced in |Me] and [Me2] . and a 
good introduction may be found in [Grj . 

The analogous structure for the short-time heat kernel is well-understood: 

Theorem 3. Fort < I, the heat kernel H M (t,x,y,x f ,y') is polyhomogeneous conormal on 
[M 2 c {x, y, x', y') x [0, 1]^; {Vi = 0, x = x' , y = y'}]. 

Moreover, there is infinite- order decay at all faces except the scattering front face sc and the 
face F obtained by the final blowup. 

Theorem [3J follows immediately from the work of Albin |Alb] ; however, the precise state- 
ment above does not appear in the literature. Therefore, in the Appendix, we give a simple 
proof using the machinery developed in |Alb] . Combining Theorem [2] with Theorem [3J gives a 
complete geometric-microlocal description of the structure of the heat kernel. This structure 
is illustrated in Figure [1] the short-time structure is the left-hand side of the diagram, and 
the long-time structure is the right-hand side. We also indicate the leading order of the heat 
kernel at each of the boundary hypersurfaces, in terms of \ft at t = 0, w = t~ 1 / 2 at t — oo, 
and x or x' respectively at all the finite-time boundaries. 

As an example of how polyhomogeneous structure may be used to read off the behavior 
in all asymptotic regimes, let a be a parameter, and fix (y,x',y'); consider the heat kernel 
H z (a 2 , a -1 , y, x', y') as a approaches infinity. In the compactified space in Figure [TJ as a 
approaches infinity, the arguments approach a point in the center of the face lbo, where the 
leading order of the heat kernel is n. Since w = a -1 , and w is a boundary defining function 
for lbo, we conclude that H z (a 2 ,a~ 1 ,y,x',y') has a polyhomogeneous asymptotic expansion 
in a~ l as a — > oo, with leading term C n a~ n . A similar analysis may be performed in any 
asymptotic regime. 

As an application, Theorems H] and [3] give us precisely the polyhomogeneous structure we 
need to define and investigate the renormalized zeta function on M: 
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Figure 1. Asymptotic structure of the heat kernel on M. 

Theorem 4. Let M be asymptotically conic. The renormalized zeta function, defined for- 
mally by $5\), is well-defined and has a meromorphic continuation to all of C. 

We may then define the renormalized determinant of the Laplacian on M by 

log( i? detA M ) = -\; f (0), 

where R (' M (Q) is the coefficient of s in the Laurent series for r (m(s) around s = 0. 

In a companion paper [S2j . we use these theorems to analyze the behavior of the deter- 
minant of the Laplacian on a family of manifolds degenerating to a manifold with conical 
singularities. We expect that this work will have applications to spectral theory and to index 
theory on singular spaces, including the study of the Cheeger-Miiller theorem on manifolds 
with conical singularities. 

1.2. Outline of the proofs. The usual geometric-microlocal approach to the fine structure 
of the heat kernel is a direct parametrix construction, which involves the construction of an 
initial approximation to the heat kernel and then the removal of the error via a Neumann 
series argument. This is the method adopted in [Alb]. However, parametrix constructions 
are not well-suited for analysis of the long-time heat kernel; the problem is global rather 
than local. In order to obtain the asymptotic structure of the heat kernel at long time, we 
instead take an indirect approach. Recall that the functional calculus shows that the heat 
kernel and the resolvent are related by 

H M (t) = Je Xt (A M + X)- 1 dX, (6) 

where T is a contour around the spectrum. In a series of papers [GHll IGH2] , Guillarmou and 
Hassell have analyzed the asymptotic structure of the resolvent (Am + A) -1 at low energy, 
again giving a complete description in all regimes. They have shown: 

Theorem 5. [GHlJ Suppose that M is an asymptotically conic manifold of dimension n > 3. 
Then the Schwarz kernel of (Am + e^k 2 )' 1 is polyhomogeneous conormal on M| sc for each 
6 G (— 7r, 7r), with a conormal singularity at the spatial diagonal and all coefficients smoothly 
depending on 9. It decays to infinite order at the faces lb, rb, and bf, with leading orders at 
sc, bf , rb , lb and zf given by 0, n — 2, n — 2, n — 2, and respectively. 
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Note that Guillarmou and Hassell require n > 3. In Section 4, we adapt the methods of 
[GH1] to extend Theorem [5] to the two-dimensional case: 

Theorem 6. Theorem^ also holds when n = 2; all the leading orders are the same, except 
that we have logarithmic growth instead of order at zf. 

In Section 2, we use geometric microlocal analysis, in particular Melrose's pushforward 
theorem, to prove Theorem [2j The key is to push the structure of Theorems [5] and [6] through 
the contour integral ([6]). The main technical theorem, which we state in slightly greater 
generality, is the following: 

Theorem 7. Let M be an asymptotically conic manifold, and let E be a vector bundle over 
M. Let A : C°°(E) — > C~°°(E) be a pseudodifferential operator with the following properties: 

a) a(A) C [0, oo]; 

b) The Schwarz kernel of the resolvent {A + e ld k 2 )~ l is polyhomogeneous conormal on Mf sc 
for each 9 G (— ir, tt), with a conormal singularity at the spatial diagonal and all coefficients 
smoothly depending on 9. Moreover, it decays to infinite order at the faces lb, rb, and bf, 
with index sets at sc, bf , rb , lb and zf given by R sc , Rbf , R r b , Rib , an d R z f respectively; 

c) The resolvent at high energy is polynomially bounded in k along with all derivatives. 
In particular, there exist ip G (ir/2,ir) and N G N with the following property: for any b- 
differential operator Q in k~ l and any scattering differential operator P on M , there is a 
constant C such that 

\PQ(A + e w k 2 )- l (x,y,x',y')\ < Ck N 

for all 9 G [—p, ifi] and all k > 1. 

Then for t greater than any fixed T > 0, the kernel of e~ tA is polyhomogeneous conormal 
on M 2 sc , where w = t~ l l 2 . It decays to infinite order at lb, rb, and bf, and has index sets 
at sc, bfo, rb , Ibo, and zf which are subsets of R sc , Rbf + 2, R r b Q + 2, Rib + 2, and R z f + 2 
respectively. 

Once we have proven Theorem [?l Theorem [2] is an almost immediate consequence, though 
there is a slight twist involving the leading orders. 

In section 3, we use Theorem [2] and some additional geometric microlocal techniques to 
analyze the renormalized heat trace and prove Theorem |H We also analyze the renormalized 
zeta function and determinant in the special case where M is exactly conic (or Euclidean) 
outside a compact set. In section 4, we extend Guillarmou and Hassell's work in |GHlj to 
prove Theorem [61 Finally, in the Appendix, we use the framework in [AlbJ to prove Theorem 

El 
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2. From resolvent to heat kernel 
The goal of this section is to prove Theorems [7| and [2j 



2.1. Preliminaries. We first give a brief summary of the key relevant concepts in geometric 
microlocal analysis; again, a self-contained introduction may be found in |Grj . A manifold 
with corners of dimension n is a topological space which is locally modeled on x M. n ~ k for 
some k; a simple example is the n-dimensional unit cube. Blow-up is a way of creating new 
manifolds with corners from old ones, and is used to resolve certain geometric singularities. 
The idea is to formally introduce polar coordinates around a submanifold of a manifold with 
corners, in order to distinguish between directions of approach to that submanifold. For 
example, consider the origin as a submanifold of To blow up the origin, we introduce 
polar coordinates (r, 6), which corresponds to replacing the point (0, 0) with a quarter-circle, 
which corresponds to the inward-pointing spherical normal bundle of (0,0) C R+. See |Mej . 
{Me2j . |Grj . and/or the appendix of [S] for a more detailed explanation and more general 
examples of blow-ups. 

By Taylor's theorem, smooth functions on a manifold with corners are precisely those 
functions which have Taylor expansions at each boundary hypersurface and joint Taylor 
expansions at every corner. Polyhomogeneous conormal distributions, which we abbreviate 
as phg or phg conormal, are a generalization of smooth functions. In particular, if we let x 
be a boundary defining function, we allow terms of the form x s (\ogx) p for any x G C and any 
p G No to appear in the asymptotic expansions at the boundary and in the joint expansions 
at the corners. The index set of a phg conormal distribution u at a particular boundary 
hypersurface H is simply the set of (s,p) which appear in the asymptotic expansion of u at 
H. 

Polyhomogeneous conormal functions are well-behaved under addition and multiplication, 
but also under more complicated operations, namely pull-back and push-forward. To discuss 
these, we first need to discuss properties of a map / : W — > Z between manifolds with corners. 
Roughly, we say that / is a b-map if it is smooth up to the boundary and product-type near 
the boundary in terms of the local coordinate models (see |Gr] for a precise definition). If 
additionally / does not map any boundary hypersurface of W into a corner of Z, and / is 
also a fibration over the interior of every boundary hypersurface, then we call / a b-fibration. 
We now have the following two theorems, due to Melrose |Me[ IMe2] , which will be critical 
in the analysis to follow: 

Proposition 8 (Melrose's Pull-Back and Push- Forward Theorems). Let f : W — > Z be a 

smooth map of manifolds with corners. 

a) If f is a b-map and u is phg conormal on Z , then f*u is phg conormal on W . Moreover, 
the index sets of f*u may be computed explicitly from those of u and the geometry of the 
map f . 

b) If f is also a b-fibration, v is phg conormal on W, and f*v is well-defined (the pushfor- 
ward is integration along the fibers, which may not converge), then f*v is phg conormal on 
Z , and again the index sets may be computed explicitly. 

Finally, we need to consider distributions which have pseudodifferential-type conormal 
singularities at submanifolds in the interior of a manifold with corners. 
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Figure 2. First blowup in creation of M| 



Definition. [Grj Let y = (xi, . . . , Xk) and z = (a^+i, . . . , x n ). A distribution u on W 1 has a 
conormal singularity at N of order m if it can be written 



with each coefficient a m ^j smooth in y and ||r . 

This definition may be extended, by using the local coordinate models, to define distri- 
butions with a conormal singularity at any p-submanifold of a manifold with corners; a 
p-submanifold is a subset which, in each local coordinate chart, may be identified with a 
coordinate submanifold. Variants of the pullback and pushforward theorems also hold for 
polyhomogeneous conormal distributions with interior conormal sing ularities [Me2l IEMM] . 

2.2. The space M| sc . We now introduce the space M| sc , which first appears in an unpub- 
lished note of Melrose and Sa Barreto [MSB], and later in (1111. To construct M| sc , we 
begin with the space Mf = [0, 1] x M x M; coordinates on this space are (k,x,y,x',y'). 
There are three boundary hypersurfaces: {k = 0}, which we call zf, {x = 0}, which we call 
lb, and {x' = 0}, which we call rb. 

First we blow up the corner {x = 0,x' = 0,k = 0}, which corresponds to the introduction 
of polar coordinates near that corner; we call the front face of this blowup bfo. We then blow 
up three codimension-2 submanifolds: we blow up {x = 0,x' = 0} and call the new face bf, 
we blow up {x = 0, k = 0} and call the resulting face lb , and we blow up {x' = 0, k = 0} 
and call the resulting face rbo- The resulting manifold with corners is shown in Figure [U 
with y and y' suppressed. Using the definition of a "b-stretched product" from the work 




where a is a classical symbol; that is, a has asymptotics as |£| — )■ oo 



3=0 
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Figure 3. The space M% . 

of Melrose- Singer [MS] , it is precisely X$(x,x', k) x N y x N y i. We use these b-stretched 
products from [MSJ throughout the arguments. 

Finally, consider the intersection of the closure of the interior diagonal {x = x' > 0, y = y'} 
with the face bf. In coordinates, this is {(x/x 1 ) — l,y = y', (x/k) = 0}, and it is marked 
with a dotted line in Figure [2J We blow this up to create a new boundary hypersurface, 
which we call sc (for "scattering"). The resulting space is M| sc , and it has eight boundary 
hypersurfaces, illustrated in Figure [3j The spatial diagonal D^ sc is defined to be the closure 
in M| sc of the interior spatial diagonal {x = x' > 0, y = y'}; its intersection with the 
boundary is marked with a dotted line in Figure |3j 

We now describe some useful coordinate systems on sc . Near the intersection of zf, 
rb , and bf , we use the coordinates 

, x f f k . 

^ X,CT = ~x ,V,y ,K = x 

In these coordinates, x is a boundary defining function (bdf) for bfo, a is a bdf for rbo, and 
k' is a bdf for zf. Similarly, near the intersection of zf, lbo, and bfo, we use the coordinates 

iii % i ™\ 
^ X ,a = ~x'' V ' V ,K = x 
Coordinates near sc are slightly more complicated; before the final blowup, good coordi- 
nates are (x',a',y,y', (x/k)). After the blowup, we use the coordinates 

v u 1 1 n v i,y~y' \ x i 

X = k(— ),Y = k , A = T ,y,k. 

x x x k 

These are valid in a neighborhood of the intersection of Dj^sc with sc and bfo; however, they 
are not good coordinates as we approach bf. 

2.3. Proof of Theorem [7]. Let A be an operator satisfying hypotheses a)-c) of Theorem [7J 
and let R(X, x, y, x', y') be the Schwarz kernel of (A + A) -1 . The spectrum of A is [0, oo], so 
i?(A, x, y, x', y') is holomorphic outside the non-positive real axis. For any a > 0, let Y a be 
the path in C consisting of two half-rays along 6 = —tp and 6 = <p, connected by the portion 
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of the circle of radius a from 9 = —<p to 9 = (p, and traversed counterclockwise. Moreover, 
let T aj i be the portion of T a along the circle of radius a, and let T a ^ be the remainder; that 
is, the two half-rays. These contours are illustrated in Figure HI 

Let F(w,x,y,x',y') be the heat kernel at time t = w~ 2 . Then by the functional calculus, 
we have 

F(w, x, y, x', y') = F(w, x, y, x', y') = / e x/vj2 R(X, x, y, x', y') d\. (7) 



r 



We let a = w 2 and T = T w 2, and then consider the integral ([7]) over T w 2 )1 and T w 2 2 separately. 
On r„,2 1, A = w 2 e ld , so dX = w 2 d9, and we have 



w 2 



r ' e £lt 'R(6,w,x,y,x',y') d9. (8) 

J — (p 

By condition b), for each 9, the integrand in (jSJ) is phg conormal on M 2 sc with a conormal 
singularity at A^ )SC , and the dependence of all coefficients on 9 e [— <p, <p] is smooth. There- 
fore, the integral (jSJ) is phg conormal on M 2 sc , with a possible conormal singularity at D WySC . 
The index sets of OH]) on M 2 sc are those of w 2 plus those of R(w, x, y, x', y'). The function w 2 
is smooth and has order 2 as a function at zf, bfo, rb , and lb and order everywhere else, 
so we add 2 to the index sets of the resolvent at those faces. This procedure gives precisely 
the index sets claimed in Theorem [71 

It remains to consider the integral over T w 2 2 - ^ w 2 ,2 consists of two half-rays; we consider 
only the half-ray corresponding to 9 = ip, as the other is analogous. Since 9 is fixed, we 
suppress 9 in the notation; the integral runs from r = w 2 to r = oo. After changing variables 
from r to s = yfr, the T w 2 2 portion of (J7J becomes: 

2 se ^ os ^ s2 / w2 e i ^ m ^ s2 / w2 R(s,x,y,x',y') ds. (9) 

First consider the behavior of the integrand as s — > oo. By condition c), R(s,x,y,x' ,y') 
is polynomially bounded together with all derivatives as s -)■ oo. Moreover, since cos (p < 
and w is bounded above, the term exp ((cos <p) s 2 /w 2 ), and hence the entire integrand of (jHJ), 
decays exponentially together with all derivatives as s — > oo. 
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Figure 5. Decomposition of R. 



In order to analyze fl5J), we break R(s,x,y,x',y') into pieces. First we separate out the 
conormal singularity in a neighborhood of D ssc and analyze using explicit local coordi- 
nates. Then we deal with the remainder, which is smooth on M^ sc , by breaking it into two 
pieces and applying the pushforward theorem. 

2.4. The conormal singularity. Using a partition of unity, we let R = R s + R C} where R c 
is supported in a neighborhood of D s sc and R s is supported away from D ssc , as in Figure 
5.4. In a neighborhood of scnbf , we use the coordinates 

(X = s( -),Y = s ,y, X = -,s). 

iAj lAj JU O 

In these coordinates, R has a conormal singularity at {X = Y = 0}. On the other hand, in a 
neighborhood of bf nzfnD SjSC , we use a slight modification of the coordinates in the previous 
subsection: 

(X = (l-^),Y = y-y',y,-,x). 
x' x 

Since we use different coordinates in different regimes, write R c = R\ + R 2 by using a smooth 
partition of unity near x/s — 1; R\ is supported away from sc (say s/x < 2) and R2 away 
from zf (say x/s < 2). This decomposition is illustrated in Figure El 
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First look at Ri. Using the explicit symbolic form of a conormal singularity, we may write 
Ri~ e^>^«^a/,x,4)|(| 2 ^(. (10) 

This is an asymptotic sum, modulo smooth functions on M S 2 SC ; we pick a particular repre- 
sentative which is supported in a small neighborhood of D S)SC , and absorb the remainder 
into R s . The coefficients cij are phg conormal in x and s/x with index sets independent of j; 
they are also smooth in y and ||r. We plug ([TO]) into ([9]) and then interchange the convergent 
s- integral with the asymptotic sum and the oscillatory integral over M. n . The result is 

Jl" ■ Jw X Isl 

By Melrose's pullback theorem |Me2j . the pullback of each aj to X^(s, x, w)xN y x S^h via 
projection is also phg conormal with index sets independent of j. As a result, the integrand 
in 



2 p i V , s £ 



e-M«^ aj {-,x } y } ±r)d8 (12) 

is phg conormal on X^(s,w,x) X N y x S^i, with a cutoff singularity at s/w = 1 and 
with infinite-order decay as s — > 00, together with all derivatives, independent of u> < 
1. Integration in s is a b-fibration from Xjf(s,w,x) to X b 2 (w,i), by Proposition 4.4 of 
[MS] , and is hence also a b-fibration when we take the direct product with iVj, x S^h. 
Moreover, integration in s is transverse to the cutoff singularity at s/w = 1. Therefore, by 
the pushforward theorem with conormal singularities (from the Appendix of |EMMj ). ( [T2[) 
is phg conormal on X%(w,x) x N y x S^l, with index sets independent of j. 

Since < s < 2x on the support of Ri, ffT^]) has a phg conormal expansion in (w/x,x). 
Therefore, the integral ( [TT[) . in the coordinates (A, Y, y, w/x, x), corresponding to the R\ 
piece of ([9]), is phg conormal in (w/x,x) with index sets independent of j, and smoothly 
dependent on y, with an interior conormal singularity at X = Y = 0. Thus ( ]TT[) is phg 
conormal on M 2 sc with a conormal singularity at the diagonal. 

We now consider R2; the analysis is similar. Write 

eWW^f^bj^yMltr'dt, (13) 

where the 6j are phg conormal in - and s with index sets independent of j, and also smooth 
in y and ||r. 

It is helpful to consider the regimes w > x/2 and w < 2x separately. First assume that 
w > x/2, and let X = w(~ — —), Y = w^ 2 ^-, and A = We expect a conormal singularity 
at X = Y = in this regime. Noting that (X,Y) = (s/w)(X ,Y), we change variables in 
( TT3[) and let C = ( S / W )C- The result is: 

„ 00 >. 

i? 2 ~ / e W>^^^6^-, S ,y )T ^)( S /^r 2 -1C| 2 -^C- (14) 

</R™ J=0 s IS. I 
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As before, plug (fT4|) into (j9]) and interchange the sums and convergent integrals: the part of 
(Q coming from R 2 is 

I e'^^fv re-^ 2 ^b^,s,y,h(s/wy-^ ds)\(\^dC (15) 
Jm n j=0 Jw s Is I 

Consider the coefficients 

<'M^ bj ^S t y C ){s/w y-2-n dg _ (16) 

If we can show that the coefficients (1T61) are phg conormal in (— ,w) with respect to some 
index sets independent of j, with smooth dependence on y and X, then ( lisp is phg conormal 
on M 2 sc with a conormal singularity at the diagonal when w > x/2. To show this phg 
conormality, note that again the integrands in ( fl6l) are each phg conormal on X^(s,x,w) x 
N y x S™7i. Moreover, the index sets are independent of j, as s/w > 1, and there is always 
infinite-order decay at s — oo, independent of w. As before, we use the pushforward theorem 
to integrate in s, and we conclude that the coefficients ( fl6l) are phg conormal on X 2 (x, w) x 
N y x with respect to index sets independent of j. Since w > x/2, this yields expansions 
in ~ and w, which is precisely what we need. 

On the other hand, suppose that w < 2x. Then we expect a conormal singularity at 
X = Y = 0. Since (X,Y) = (s/x)(X,Y), we change variables in ( fl3l) . letting (' = (s/x)£. 
Following the exact same procedure as in the w > x/2 case, we see that the part of (jUj) 
coming from when w < 2x, is: 

f e i(x,ry(CiA)J2( fe-W^6/, S) 4)(^ ^IC'I^VC- (17) 

The coefficients can be rewritten as 

{w/xy~ 2 ~ n j°° e-M^bjQ s , y, l)( s / w )^- d s . (18) 

But these are just {w / xy~ 2 ~ n times ([TBI) . (w/xy~ 2 ~ n is phg conormal on A^(w, x), so (fT8|) is 
also phg conormal on X%(w, x) for each j. Since we are only considering w < 2x, the orders 
only improve as j increases. In particular, all the coefficients are phg conormal on X 2 (w, x) 
with respect to subsets of the index set of the j = coefficient. This is again sufficient to 
prove that (TTT]) . and hence the part of (J9j) corresponding to R c , is phg conormal on M 2 sc 
with a conormal singularity at the diagonal. 

Technically, we need to compute the index sets of the coefficients of the conormal singular- 
ity at each boundary face of M 2 sc . This may be done directly via the pushforward theorem, 
but it is easier to apply the analysis we will develop in the next section. Note that aj and bj 
may be viewed as phg conormal functions on the diagonal D s sc C M 2 SC , with fixed index sets 
R sc , Rbf , and R z f at the boundary hypersurfaces sc, bfo, and zf. Observe that they can be 
extended smoothly to functions defined in a neighborhood of D s sc which are themselves phg 
conormal on M 2 sc with the given index sets; call these extensions Cj. Then the coefficients 
of the conormal singularity of the integral (Q are the restrictions to the diagonal of 

oo 

2se- (s/w)2e ^c 3 (s,x,y,x\y') ds. 
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Applying Lemma [9] to each Cj, these coefficients are all phg conormal on M 2 sc , with index 
sets obtained by adding 2 at the faces bfo, rbo, lbo, and zf. Therefore the restrictions to the 
diagonal are all phg conormal on D WySC , with leading orders matching those in Theorem [7J 
as expected. This completes the analysis of the conormal singularity. 

2.5. Finishing the proof. It remains to consider the integral 

POO 

/ 2se-( s / w ^R s (s,x,y,x',y') ds, (19) 

J w 

where R s (s,x,y,x',y') is phg conormal on M 2 SC and smooth across the diagonal. We claim: 

Lemma 9. Let T(s,x,y,x',y') be any function which is phg conormal on M 2 sc , smooth in 
the interior, decaying to infinite order at lb, rb, and bf, and has at most polynomial growth 
in s as s — )■ oo along with all derivatives (as in condition c) of Theorem^). Then 

POO 

/ 2se^ s/w)2eliP T(s,x,y,x',y') ds (20) 

J w 

is phg conormal on M 2 sc . Moreover, if the index sets of T at the various boundary hyper- 
surfaces are T sc , T b f , T z j, T rbo , and Ti bo , then the index sets of [Wty are: 

T sc , T b f + 2, T z f + 2, T r b + 2, T^ + 2. 

We defer the proof for the moment. Applying Lemma to T = R s , we conclude that (1191) 
is phg conormal on M 2 sc with index sets precisely as in Theorem [7J Combining this with our 
analysis of R c , we have now shown that F(w,x,y,x',y') is phg conormal on M 2 SC possibly 
with a conormal singularity at the spatial diagonal, and with leading orders as specified in 
Theorem [7J However, F(w,x,y,x',y') is a heat kernel, so it has no conormal singularity at 
the diagonal. This completes the proof of Theorem [7J 

Finally, to prove Theorem [21 we apply Theorem [7J Condition a) is true since the Laplacian 
is essentially self-adjoint and non-negative. Condition b) follows from Theorems [5] and [6j 
Finally, condition c) is a well-known consequence of the semiclassical scattering calculus. The 
scattering calculus was first introduced by Melrose in |Me3] and the semiclassical version was 
developed by Vasy, Wunsch, and Zworski among others |VZ| IWZ] . The exact statement we 
need, along with a summary of the semiclassical scattering calculus, may be found in section 
10 of [HWj . Applying Theorem [7J gives us the polyhomogeneity we claim, and once we plug 
in the leading orders from [GHlJ and the Appendix, we see that the heat kernel has leading 
orders of at sc and n at each of bf , rb , and lb . 

Unfortunately, Theorem [7J does not by itself give us the claimed order-n behavior at zf; 
instead, we only see quadratic decay at zf when n > 3 and decay of the form w 2 log w when 
n = 2. However, by the estimate of Cheng-Li- Yau (Theorem [1]), the heat kernel is uniformly 
bounded for large time by Ct~ n l 2 = Cw". Thus the leading order of the heat kernel at zf 
must actually be at least n, which completes the proof of Theorem [2J Notice that the order 
calculation in Theorem [7J is not sharp in this case. We expect that the lack of sharpness 
is due to some complex-analytic structure of the resolvent which is lost in our real-analytic 
approach. 

2.6. Proof of Lemma [91 We now prove Lemma [9j the proof involves extensive use of 
Melrose's pullback and pushforward theorems. First, write T(s, x, y, x', y') as 7\ + T 2 , where 
Ti is supported away from sc and T 2 is supported in a neighborhood of sc. This partition is 
illustrated in Figure El Then decompose (|2"0"]) into two integrals, corresponding to Ti and T 2 . 
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: T* 




Figure 6. Decomposition of T. 

Consider the first integral: 

( r X{{s > w^wxx'e-W^T^x^x^y')-) ( 2 1) 
\J s J w x x' 

Notice that Ti(s,x,y,x',y') is phg conormal on M 2 sc but supported away from sc, so it 
is in fact phg conormal on the blown-down space X^(s,x,x') x N y x N y ,. The rest of 
the terms in the integrand are phg conormal on X 2 (s,w), with a conormal singularity at 
s = w. By the pullback theorem, we conclude that the integrand is phg conormal on 
X£(s,w,x,x') x N y x N y r, with a conormal singularity at s = w. It decays to infinite order 
as s — > oo uniformly with all derivatives. 

Now consider the map to X^(w,x,x') x N y x N y , given by integration in s. This map is 
a b-fibration, with fibers that are transverse to the p-submanifold {s = w}, where there is a 
conormal singularity So by the pushforward theorem , (I2TI1 is phg conormal on Xjf(w, x,x')x 
N y x N y i. But X'^(w, x, x') x N y x Ny, is a blow-down of M 2 W sc (w, x, y, x' , y'), so (f2T1) is phg 
conormal on M 2 sc (w, x, y, x', y'), as desired. 

We now consider the integral 

( f°x({s > W })2 S We-^) 2 ^)T 2 ( S ,x, 1/) x', 2/ ')-) — (22) 
J s w x x' 

Let a = (f? — 1, y — y); a is an n- dimensional coordinate, and M 2 S sc is created from 
Xfr(s,x,x') x N y x N y , by blowing up {a = ^ = 0}. In particular, T 2 (s,x,y,x',y'), having 
compact support in x/x', is phg conormal on 

[X 2 b (s,x)xN yi {a=^ = 0}]. 

This space is the subset of M 2 sc with 1/2 < x/x' < 2, so label its boundary hypersurfaces 
bf, sc, bfo, and zf. In this labeling, T 2 is supported away from zf, decays to infinite order at 
bf, and has leading orders t sc at sc and tbf at bf . 
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Figure 7. The space S. 

We analyze the integrand in ( 122|) as a function on the space 

S = (X^(s,w,x) Hs>w)x B(a) x N y . 

Here B is the unit ball in W 1 . A diagram of 5* is given in Figure [7J with a and y suppressed; 
we label the boundary hypersurfaces A-E. 

We now define an iterated blow-up of S. Let Pi be the p-submanifold of S given by 
An{<7 = 0}. Blowing up Pi creates a new space Si = [S; Pi]; call the front face of this blowup 
F. Now let P2 be the p-submanifold of Si given by the closure of the lift of D° fl {a = 0}. 
Then let 

S 2 = [Si;P 2 ] = [[S;Pi];P 2 ], 

and let G be the new front face. The following two propositions allow us to analyze (|22|) ; 
their proofs are deferred for the moment. 

Proposition 10. The map 

n w :S 2 ^ [(X^s,x) n s > x) x N y ; {a = ^ = 0}], 

given in the interior of S 2 by projection off the variable w and extending continuously to the 
boundary, is a b-map. 

Proposition 11. The map 

7r s :Si^ [Xl( W) x) x Bi(a x N y , {a = ^- = 0}], 

given in the interior of Si by projection off the variable s and extending continuously to the 
boundary, is a b-fibration. Moreover, if we let pn be a bdf for each hypersurface H, we have 

M*(P*f) = PCPE, (7T s )*(P6/ ) = pBpD, ^sT(psc) = pF, {^s)*{pbf) = P A- (23) 
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Since T 2 is supported in {s > x}, Proposition [10] and the pullback theorem imply that the 
pullback of T 2 is phg conormal on S^. Moreover, the remainder of the integrand in fj22|) is phg 
conormal on (X%(s, w) fl {s > w}), so pulling back first to X^(s, w, x) fl {s > w}, and then to 
S2, we see that it is phg conormal on S2 as well. Therefore, the entire integrand in (I2"2"j) is phg 
conormal on However, since it is phg conormal on S2 = [Si; P2] and vanishes to infinite 
order at the front face G, it is actually phg conormal on Si. By the pushforward theorem 
from [EMM], since ir s is a b-fibration transverse to the conormal singularity at s = w, the 
pushforward ff22l) is phg conormal on the target space [X$(w, x) x Bi(a x N y ; {a — ^ = 0}]. 
From Figure [9] we see that this space is a subset of sc ; we have therefore shown that 
( 122|) is phg conormal on sc . This completes the proof of the polyhomogeneity statement 
in Lemma [9] modulo the proofs of Propositions [10] and [11] 

It remains to check the index sets claimed in Lemma [9] However, this calculation is a 
straightforward application of the pullback and pushforward theorems (explicit descriptions 
of the pullback and pushforward index sets may be found in |Grj ). A computation of the 
leading orders may be found in [S] and computing the index sets themselves is no harder. 

2.7. Propositions [TU1 and [TTl Finally, we prove Propositions [TU1 and ITT1 The proofs involve 
explicit local coordinates. For simplicity, redefine 5* to be Xjj(s,w,x) fl {s > w} x B(a); 
there is an extra N y , but it is not involved in any of the blow-ups. Divide S into regions 
I- 1 1 1 , where s > w > x in region I, s > x > w in region II, and x > s > w in region III. This 
partition is illustrated in Figure [8] Projective coordinate systems in regions I-III are given 
by 

(/):(*,-,-, a), (II) :(*,-,-,*), (III) : (x, -, -, a). (24) 

S W S X X s 

Let pa, ■ ■ . , Pe represent boundary defining functions for the boundary hypersurfaces. Each 
coordinate corresponds to a local bdf in the specified region, or to a; for example, in (II), 
s = p B ,- s = Pd, and f = p c . 

The new space Si is created by blowing up An{a = 0}. This blowup is localized to region 
I, so pa — — ■ Call the new face F. To find coordinates on Si, we split region I into two new 
regions: IV, where \a\ < ~, and V, where — < \a\. Then the coordinates (I) are replaced 
with two new coordinate systems: 

(IV) : (s,-, -,——), (V) : [s, \a\). (25) 
s w [x/w) s \a\ \a\ 

Proof of Proposition Q7] 

We now prove Proposition [TT] We must show that projection off s, extended to the 
boundary in Si, is a b-fibration onto [X%(w,x); {a = — = 0}]. This space is illustrated in 
Figure O the faces are labeled sc, bf, bf , and zf, to illustrate that the space is a subset of 
sc . We split the space into three regions along the dotted lines in Figure M, Region 1 is 
a neighborhood of zf, Region 2 is near the interior of sc, and Region 3 is the remainder (a 
neighborhood of bf). Coordinates on [X$(w,x); {a = — = 0}] are: 

. . , w . /N / x a . /„s / (x/w) a . . 

1 :x,-,a, (2):(w,- TT -A, 3 : (w, ±LJ-, -- ,<t . 26 
x w (x/w) \a\ \a\ 

We first check that tt s is a b-map. This can be done using the local coordinate models. 
For example, region 1 is projected from regions II and III. Examining the map from region 
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Figure 8. Decomposition of S. 
i 




Figure 9. [X*(w,x); {a = * = 0}]. 



II to region 1, we see that 



W X w 

{x,—,a) = {s- -,— ,<t). 

X S X 
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Hence every coordinate in region 1 is a product of (positive integer powers of) coordinates 
in region II. The same is true for region III, as 

( w -\ ( w s -\ 

[ x , — ,cr) = (x, ,cr). 

X S X 

Therefore ir s is a b-map onto region 1. Using this type of analysis for regions 2 and 3, it is 
easy to see that ti s is a b-map. 

To check that n s is a b-fibration, we look at the pullbacks of the boundary defining 
functions by using the local coordinate models. For example, in region 1, p z f = — . The 
pullback to region II is — = pc, and the pullback to region III is — • - = PcPe- These are 
consistent, as region II is bounded away from E. So the pullback (^ s )*(Pzf) = PcPe- We can 
perform this analysis on each face of the target space, and obtain: 

M*(Pzf) = PCPE, (7T S )*(P6/o) = PBpD, (^s)*(psc) = pF, (^s)*(pbf) = pA- (27) 

We conclude that 

7T S (C) = 7C S (E) = Zf, 7C S (B) = 7T S (D) = bf , n s (F) = SC, 7f s (A) = 

In particular, no boundary hypersurface is mapped into a corner, so the first condition in 
the definition of b-fibration holds. As for the second (fibration) condition, the map from 5* 
to X^(w, x) x B(a) given by integration in s is a b-fibration, by Propositions A. 10 and A. 12. 
Since the blow-up only involves faces A and F, we conclude that the fibration condition is 
satisfied on the interiors of sc and bfo- 

Finally, near A and F, we have coordinates (IV) and (V). The map from (IV) to (2) 
fixes the last two coordinates and maps (s, — ) to w = s ■ —. Since the map from X^(s,w) 
to R+ is a b-fibration, so is the map from (IV) to (2). Analysis of the map from (V) to (3) 
is identical, so n s is a b-fibration near A and F, and hence everywhere. This completes the 
proof of Proposition [TT1 

Proof of Proposition [70 

To prove Proposition [1XJ, we need to make another blowup, of the submanifold P<i given 
by the closure in S\ of D° PI {a = 0}. As before, call the new face G. This new face intersects 
regions II and IV and requires new coordinates in each. In region II, we create G by blowing 
up = 0, a = 0}, and in region IV, we create G by blowing up = 0, = 0}. 

Therefore, we split region IV into regions VI and VII, where jfp^ < 7 and j^p^ > 7 

respectively; note that = 7^y- Also split region II into regions VIII and IX, where 

a < - and - < a respectively. 

Then there are six different local coordinate systems on S2: (HI), (V), and: 

(VI) ■ (s, -, -, -— - , (VII) : (s, -, i—r^, tzt), 

s w (x/s) \a\ w [x/w) \a\ 

i\r j j j\ f x w a (x/s) w a 

(VIII) : (s, -, -, (IX) : (s, —— , -, a ). 

s x (x/s) \a\ x \a\ 

We want to show that ti w is a b-map onto [(Xjj(s, x) D s > x); {a — - — 0}]. Since s > x, 
we can ignore region III. Following our analysis of the target space in Proposition [TTJ we 
break [(X$(s,x) fl s > x); {a = - = 0}] into two regions, region 2 (where \a\ < -) and 
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region 3 (everywhere else). Coordinates on these regions follow the previous analysis, with 
w replaced by s: 

(2):(.,f,^j), (3):( S ,<^,^). (28) 

At this point, it is simply a matter of checking the local coordinate models to prove that 
tt w is a b-map. Note that regions VI and VIII map to region 2, whereas regions V, VII, and 
IX map to region 3. The details may be easily written down using the local coordinates as 
in the proof of Proposition [TTJ 

It is worth mentioning that tt s is not a b-fibration from S2 to [X%(w, x); {a = — = 0}]; G 
gets mapped into the intersection of bfo and a — 0. 

3. Renormalized heat trace and ZETA FUNCTION 

In this section, we define the renormalized heat trace, zeta function, and determinant 
on an asymptotically conic manifold M. The first step is to define the renormalized trace. 
This definition is inspired by Melrose's b-heat trace, which is a renormalized heat trace for 
manifolds with asymptotically cylindrical ends. Albin also defined renormalized heat traces 
in the asymptotically hyperbolic setting [Alb] ; later, Albin, Aldana, and Rochon defined 
and investigated a renormalized determinant of the Laplacian on asymptotically hyperbolic 
surfaces |AARj . 

3.1. The renormalized heat trace. Pick any cutoff function Xi( r ) 011 ^+ which is sup- 
ported on {r < 2} and equal to 1 on {r < 1/2}. Assume that either 

a) Xi( r ) is a non-increasing smooth function of r (smooth cutoff), or 

b) Xi( r ) is precisely the characteristic function of [0, 1] (sharp cutoff). 

Then for any 5 < 1/2, let Xi,s be a function on M, equal to Xi{ r $) f° r r >1 and equal to 
1 inside {r = 1}. Consider the integral 

/ X iMH M {t,z,z) alz. (29) 

JM 

We will prove the following theorem: 

Theorem 12. Let Xi be either the smooth or the sharp cutoff. The integral ( f^Pj) has a 
polyhomogeneous expansion in 5 for each fixed t. Moreover, the finite part at 5 = 0, which 
we denote P(t), has polyhomogeneous expansions in t at t = and t -1 at t = 00. 

This theorem allows us to define the renormalized heat trace on an asymptotically conic 
manifold. Roughly, this corresponds to integrating the heat kernel on the diagonal over 
regions where r < and renormalizing by subtracting the divergent parts at 5 = 0. 
Renormalization in this fashion is often called Hadamard renormalization (for details, see 
[AibllAAR] ). 

Definition. Let Xi{ r ) be the sharp cutoff. The renormalized heat trace, denoted R TrH M (t), 
is the finite part at 5 = of ( QS| ). 

We now prove Theorem [121 The key ingredient is the following observation on the structure 
of the heat kernel on the diagonal, which is a consequence of the structure theorem we have 
proven for the heat kernel. The asymptotic structure of H M (t,x,y,x,y) reflected in this 
proposition is illustrated in Figure (TDJ 
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^X 



X 



A 



Vi= 



^ w ^ ^ ^ blow-up of {w = x = 0} 



x = 



Figure 10. Asymptotic structure of H M (t,x,y,x,y). 



Proposition 13. a) Fort bounded above, H M (t,x,y,x,y) is phg conormal in (y/i,x), with 
smooth dependence on y. 

b) For t bounded below, let w = t~ 1//2 ; then H M (t, x, y, x, y) is phg conormal as a function 
of w and x on X%(w,x), again with smooth dependence on y. 

Proof. Proposition [T3] is an immediate consequence of restricting to the spatial diagonal D 
in Theorem [2j since D is a p-submanifold of the space on which H is polyhomogeneous, 
the restriction of H M to D is also polyhomogeneous. Comparing D with Figure UHl we see 
that D is precisely the space described in Proposition [13] □ 

To prove Theorem [121 we analyze ( 129]) . which may be rewritten as 

Xi(S/x)H M (t,x,y,x,y)x' n - 1 dx dy+ H M {t,z,z) dz. (30) 

J\z\<l 

First analyze the second term; the region {\z\ < 1} is bounded away from spatial infinity. 
Therefore, by Proposition [131 H (t, z, z) has polyhomogeneous expansions in t at t = and 
t -1 / 2 , hence t -1 , at t = oo, and these expansions are uniform in z with smooth coefficients. 
Integrating in z results in a function of t which is phg conormal at t — and t = oo; this 
function contributes only to the finite part P(t) at 8 = and satisfies the polyhomogeneity 
claimed in Theorem [T2J 

It remains to analyze the first term in f l3Up . We consider the small-t and large-t regimes 
separately, analyzing the integrand 

X i(S/x)H M (t,x,y,x,y)x- n - 1 (31) 

in each regime as a function of (t, x, 5). In each case, Xi{$/ X ) is phg conormal on Xjj(x, 5); if 
Xi is the sharp cutoff, there is also a cutoff singularity, which is a type of conormal singularity, 
at S/x = 1. 

For small t, H M (t,x,y,x,y) is phg conormal in (\/t,x), so (J3TJ) is phg conormal on 
R+(Vc) x X 2 (x,<5), possibly with a conormal singularity at 8/x = 1. The projection map 
ii x is a b-fibration from this space onto the first quadrant in (yt, 8) and is transverse to 
8/x = 1; moreover, the integral in x is well-defined, as the integrand is supported away from 
the x = 0,8 > face. By the pushforward theorem from [ EMM ], the first term of (130]) is 
phg conormal in {y/t, 8) for bounded t. 

On the other hand, for large t, H M is phg conormal on X£(w,x) and Xi is P n g conormal 
on Xl(x,8). Since the maps from X^(w,x,8) to each of these spaces are b-maps (also b- 
fibrations), the integrand is phg conormal on X% (w, x, 8) by the pullback theorem; there may 
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Vt = o 



w = t~ l/2 

^ ( blow-up of {w = 8 = 0} 



5 = 



Figure 11. Asymptotic structure of (1301) . 



again be a conormal singularity at 5/x = 1. Integration in x is pushforward by a b-fibration 
onto X%(w, 5). Again, the integrand is supported in {x > 5}, and the fibration is transverse 
to 5/x = 1, so we apply the pushforward theorem from [EMMJ to conclude that the first 
term in (13T)|) is phg conormal on X%(w, 5) for bounded w. Combining these results, we have 
shown that (13"U|) is phg conormal on the space in Figure [TTJ 

In particular, for any fixed t, (I30I) has a polyhomogeneous expansion as 8 — > 0. Moreover, 
P(t) is simply the coefficient of the t° term at the (0 < t < oo, 5 = 0) face. By the definition 
of phg conormality (also see the discussion surrounding Lemma A. 4 in |Maj ) . P(t) therefore 
has polyhomogeneous conormal expansions at t = and t = oo. This completes the proof 
of Theorem I 



Using Theorem [121 we now define the meromorphic continuation of the renormalized zeta 
function: 

r (m(s) = =pr / (*TrH ii )(t)ir 1 dt. (32) 
r (s) Jo 

We break up the integral (|32|) at t = 1, and consider first the short-time piece: 

_L f\*TrH M {t))t s - 1 dt. 
r (s) Jo 

By the phg conormality of R TrH (t), we can write for any N > 0: 



(33) 



k N 



R TrH M {t) = ^a^(logt) P! + 0{t N ). 



i=0 



Plug this expansion into ( l33l) . The 0(t ) contribution is well-defined and meromorphic 
whenever 9fts > — N, and the continuations of the other terms are integrals of the form 



r s) 



t 



Zi+S— 1 



Oogt)«. 



These integrals may be evaluated directly, and give explicit meromorphic functions of s, each 
with finitely many poles. Therefore, (1331 . though initially defined only when 9ffe > — zq, has 
a meromorphic extension to all of C. 

On the other hand, the long-time piece is 



r(a) A 



TrH M (t))t s ~ l dt 



(34) 
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Writing u = 1/t and substituting, this becomes 

1 

(*T r #^(i/ w )) w (--)-i du. 

We have a phg conormal expansion for R TrH M (1/u) as it — >• 0; say the leading order term 
is of the form u 2 °° (log m) p . Proceeding exactly as in the analysis of fl33|) . we conclude that 
( l34"j) . though initially defined only when 9ft(— s) > — Zoo, has a meromorphic continuation to 
all of C. 

This allows us to define the renormalized zeta function and determinant on any asymp- 
totically conic manifold M. 

Definition. The renormalized zeta function on M, r ( ) m{s), is given by the meromorphic 
continuation of ([Hp. 

Depending on the orders, there may be no s in C for which fl32l) is defined; however, 
once we split the integral at t — 1, both pieces are defined in half-planes and continue 
meromorphically to all of C. 

Definition. The renormalized determinant of the Laplacian on M is e~ ^'m (°\ where R (' M (0) 
is the coefficient of s in the Laurent series for Cm(s) at s = 0. 

3.2. Manifolds conic near infinity. We now specialize to the case of manifolds which 
are precisely conic outside a compact set. In particular, let Z be any asymptotically conic 
manifold without boundary which is isometric to a cone outside a compact set. Without 
loss of generality, assume that Z is isometric to a cone when r > 1/2. We examine the 
asymptotic expansion of f z xi,s(z)H z (t, z, z) dz as 5 — > 0; the finite part is precisely the 
renormalized heat trace. However, for applications, such as the study of conic degeneration 
in |S2j . we are also interested in identifying the divergent terms in the expansion. The fact 
that Z is conic near infinity allows us to identify those terms: 

Theorem 14. Let Z be conic near infinity as above, and let xi be the sharp cutoff. Then 
f z Xi,s( z )H z (t, z, z) dz has the following asymptotic expansion as 5 — >■ 0: 

» n— 1 

/ X i,s(z)H z (t,z,z) dz = J2fk(t)S k ~ n + fio 9 (t)\og6 + R TrH z (t) + R(6,t). (35) 

^ Z k=0 

Here R(S,t) goes to zero as 5 goes to zero for each fixed t. Moreover, if we let Uk(l,y) be 
the coefficient oft^k — n)/2 in the short-time heat expansion on at the point (l,y), then 
fk(t) = (f N u k (l,y) dy) tlk ~yj 2 and f\ og {t) = - J N u n (l, y) dy. 

We now prove Theorem [T4l Note first that f z xi,s{z)H z (t,z,z) does in fact have a poly- 
homogeneous expansion in 5, by Theorem [T2], so it is just a matter of identifying the terms. 
The proof involves a comparison of the heat kernels on Z and on CV; Z and are identical 
near infinity, which allows us to formulate and prove the following lemma: 

Lemma 15. Let Cn be the infinite cone over n. Then \H Cn (t, z, z) — H (t, z, z)\, defined 
whenever r > 1, decays to infinite order in \z\ as \z\ goes to infinity. 

Proof. Let Z = Z fl {r > 1}. It is a complete manifold with boundary at r = 1, and is 
a subset of both Cn and Z. On Z, H N (t,z,z') — H z (t,z,z') is a solution of the heat 
equation for each z', with initial data equal to zero and boundary data at r = 1 given by 



n-(-s)) Jo 



THE HEAT KERNEL ON AN ASYMPTOTICALLY CONIC MANIFOLD 



23 



H° N (t, 1, y, z') - H z (t, 1, y, z'). Fix any T > 0. We claim that for all i < T, y G N, and z' 
with \z'\ > 2, there is a constant K so that the absolute value of the boundary data is less 
than K. 

We show this for H CN (t, 1, y, z') and for H z (t, 1, y, z') separately. For CV, by scaling and 
noting that \z'\ < 2, 

H c "(t, 1, y, z') = \z'\- n H c »(t/\zf, l/\z'\, y, z'/\z'\) < H c »(t/\zf, l/\z'\, y, 1, y'). 

For each fixed y', the heat kernel with point source at {l,y') is continuous for r < 1/2 (i.e. 
the tip of the cone), and hence is bounded for r < 1/2 and for t/\z'\ 2 < t < T by some 
universal constant K. Since y' varies only over a compact set, the proof is complete. As for 
H z (t,l , y, z'), consider the region 

W = {(t,l,y,x',y')\ t < T,x' < i y G N,y' G N} 

as a subset of (t, x, y, x', y') space. The kernel H z has infinite-order decay at each boundary 
hypersurface of the space in Figured] with which W has nontrivial intersection. We conclude 
that H z is bounded on W, so there is a constant K so that H z (t, 1, y, z') < K for all t < T, 
all y G N, and all z' with \z'\ > 2. 

Since we have an upper bound for the boundary data, we can construct a supersolution and 
apply the parabolic maximum principle. Let g(t, r) be the solution of the heat equation on Z 
with zero initial condition and boundary data at r = 1 equal to K for all t. By the maximum 
principle, we see that, uniformly for \z'\ > 2 and t < T, \H Cn (t, z, z' ,t) — H z (t, z, z')\ < 
g(t, \z'\). We claim that g{t,r) decays to infinite order in r, uniformly in t for t < T. We 
can see this either directly from Bessel function expansions or by constructing a further 
supersolution g(t,r) modeled on the heat kernel on M. n . For completeness, we write down 
g(t,r): 

, Kla JK 1 _^L_ 
9it ' r) = (4^7* 2- ( t -A:)"/2 e 4(t ~ k) Mt>k}, 

where a = (87r) n,/2 e 1//4 and T is the greatest integer less than or equal to T. This super- 
solution has the uniform exponential decay property we want, so a final application of the 
parabolic maximum principle finishes the proof of the lemma. □ 

The following corollary is an immediate consequence: 

Corollary 16. For any fixed t and any Xi,s (either a sharp cutoff or a smooth cutoff), 

I f XiAz)H° N (t,z,z) dz- f X iAz)H Z (t,z,z) dz\ (36) 

converges as 5 — >■ 0. 

It now suffices to show that 

f Xi,s( z ) HCN (t,z,z) dz 

J Cn 

has a divergent asymptotic expansion of the form claimed in Theorem [T3J as (136]) converges 
as 5 — > and hence contributes only to the finite part of the expansion. (Recall that x is 
the sharp cutoff). 
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Lemma 17. Fix t. The divergent terms in the expansion of f\ z \ K1 i S H CN (t, z, z) as 5 — > 
are given by: 

n ~ X » f(k-n)/2 



V(/ u k (l,y)dy)- 5 k - n -C\og5, 

UJn k-n 



where C is equal to f N u n (l,y) dy. 

Proof. The integral is, modulo a term independent of 5, 

/ / H Cn (t, r, y, r, y)r n ~ x dy dr. 
Ji Jn 

By the conformal homogeneity of Cn, H Cn (t,r,y,r,y) = r~ n H CN (t/r 2 , 1, y, 1, y). So the 
integral becomes: 

^ y ^(1,1,2/, l,^- 1 dydr. 
Now let s = t/r 2 and switch to an integral in s; we get: 

\ [ (/ H cN {s,l,y,l,y)s- l dy)ds. (37) 
From short-time heat asymptotics, we know that 

(/ H cN (s,l,y,l,y)dy) = J2( Mh y) dy)s^' 2 + / u n (l, y) + ©(s 1 / 2 ). (38) 
Jn k=0 Jn Jn 



Here Uk(l,y) are the heat coefficients on the cone at the point (l,y). We plug into 
(137]) and get 

"Zi /• +(fc~n)/2 /• 

£(/ u k (l,y)dy)- 6 k - n -( u n (l,y)dy)\og5 + g(5,t), 

where g(S,t) is finite as 5 — > 0. This is what we wanted to prove. □ 

Combining this lemma with the preceding corollary and the definition of the renormalized 
heat trace completes the proof of Theorem [TH 

Finally, it is also useful to investigate the analogous divergent expansion when a smooth 
cutoff, rather than a sharp cutoff, is used. 

Lemma 18. Let Xi( r ) be as in a): smooth and non-increasing, supported in r < 2 and 1 
when r < 1/2. Then: 

» 71—1 

/ X iAz)H Z (t,z,z) dz = J2 l kfk(t)S k - n + f log (t)\og6 

J Z k=0 

+ ( R TrH z (t) + l l09 f l09 {t)) + R(5,t), (39) 

where h = — / x / 2 x' 1 (r)r fc_n dr, li og = — jy 2 x[( r ) l°g r dr, and R(8, t) goes to zero as 5 goes 
to zero for every fixed t. 
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Proof. Let £(r) be any function which is equal to a constant a for r < 1/2 and supported in 
{r < 2}. For any 5 < 1/2, we may define a function £s( z ) on Z by letting £5(2) be equal to 
£<5r for r = |z| > 1 and a for {r < 1}. Then consider the integral 

/ Uz)H Z (t, Zl z) dz, (40) 
J z 

and examine its behavior as 5 — > 0. 

When £(r) is the characteristic function of [0, 1], we have the expansion (135]) . By replacing 
5 with S/b for any 6 e [1/2,2], we can compute the 5 — >• expansion of (140]) for £(r) equal 
to the characteristic function of [0,6]. By linearity, we see that the expansion of (j4T)]) for 

f (r) = (Ah)x[ a ,b] is: 

n-1 

J](A/ i )/ fc (t)(6 fe ^ - a k - n )5 k - n + (A/i)/ w (t)(log& - log a) 

k=0 

+(Ah)(R{8/b,t) - R{5/a,t)). 

Now let £(r) = Xi( r ) ~~ x{ r — 1}; this is the difference between the sharp and smooth 
cutoffs. Since the expansions of (j40"l) are linear in £, we can approximate by step functions 
and then integrate by summing over thin horizontal rectangles. Assume for simplicity that 
Xi{ r ) — 1 for all r < 1 (in the general case, there are some negative signs, but we get the 
same answer). The thickness of the rectangle at height h is Ah. The length of the rectangle 
is Xi 1 ^) — 1. Putting all of this together, the expansion of (140!) with respect to £(r) is: 

/ {{Xi\h)) k - n -l)dh)5 k ~ n + f log {t){ (\og X ^(h))dh + R(5,t), 

where R(5, t) is the contribution from the remainder terms. 

Finally, perform the change of variables u = x^ 1 (/i), then add the expansion for x{ r < 1}; 
we obtain precisely the expansion claimed in the statement of the lemma. This finishes the 
proof, as long as we can control the remainder term R(d~,t). Indeed, for each fixed t, we 
claim that R(S, t) goes to zero as 5 goes to zero; define a new function S(S, t) by letting 

S(8,t)= sup \R(6/j,t)\. 

l/2<7<2 

When £(r) = (Ah)[a,b], the remainder is bounded in absolute value by (Ah)S(S, t). So the 
integral from h = to 1 is bounded by S(S,t), which goes to zero as S — > 0; this shows 
boundedness of the remainder term and finishes the proof of the lemma. □ 

It is worth examining the dependence on the zeta function on the choice of cutoff Xi] we 
used a sharp cutoff to define it, but we could use a smooth cutoff instead. In this case, the 
finite part of the divergent (^-expansion changes from R TrH z (t) to R TrH z (t) + h gfiog(t). 
But fi g and k og are constants. So the renormalized heat trace only depends on the choice 
of cutoff function by the addition of a constant, independent of t. However, it can be 
easily shown by breaking up the integral at t — 1 that the meromorphic continuation of 
J °° Ct h t s ~ l dt is identically zero for any constant C. We have shown: 

Proposition 19. Let Z be conic near infinity. The renormalized zeta function and deter- 
minant of the Laplacian on Z are independent of the choice of cutoff function Xi,s- 
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We have now shown the existence of a renormalized zeta function and determinant of the 
Laplacian on any asymptotically conic manifold M; moreover, when M is conic near infinity, 
we have computed the divergent terms in the expansion which leads to those renormaliza- 
tions. 



4. The low-energy resolvent in two dimensions 

In this section, we extend the techniques used by Guillarmou and Hassell in [GH1] to 
prove Theorem [6j In particular, we construct the low-energy resolvent on an asymptotically 
conic surface. The resolvent is 

R(9,k,x,y,x',y') = (A M + e l9 k 2 y l (x, y, x', y'). 

For simplicity, we set 9 = 0, so that R is a function of (k, x, y, x', y'). At the end of the 
section, we return to discuss allowing arbitrary 9 G [—(p, <p] and showing smoothness in 9; 
however, this is not difficult. 

4.1. Strategy. Our goal is to construct the Schwarz kernel of the resolvent, R(k, x, y, x', y') 7 
as a distribution on M% sc . To do this, as in [GHlJ, we will first construct a parametrix G(k) 
so that (Am + k 2 )G(k) — Id+ E(k), where E(k) is an error term. G(k) will be a family (in 
k) of pseudodifferential operators on M whose Schwarz kernel is polyhomogeneous conormal 
on M 2 sc with an interior conormal singularity at the spatial diagonal. By examining the 
leading order behavior of the equation (A M + k 2 )G(k) = Id at each boundary hypersurface of 
M 2 sc , we obtain a model problem at each hypersurface. The leading order of the parametrix 
G{k) at each hypersurface should solve the model problem. We first choose solutions of the 
model problem at each hypersurface, and then check that they are consistent; that is, that 
they may be glued together to obtain a parametrix G(k). Finally, we analyze the error E{k) 
and show that it can be removed via a Neumann series argument. 

In order to define the appropriate space of pseudodifferential operators, we use certain 
density conventions, all the same as in |GH1] . We consider P = Am as an operator on 
scattering half- densities by writing 

PUix^x-^dxdy^ 2 ) = (AMf^x^lx-^dxdyl 1 / 2 . 

As in |MSB] and |GH1] . we expect a transition between scattering behavior for k > and 
b-behavior at k — 0, which leads us to define the conformally related b- metric g b = x 2 g. 
The space (M, g b ) is asymptotically cylindrical. We then define P b = x~ n / 2 ~ 1 Px n / 2 ~ 1 with 
respect to scattering half-densities. However, we want to consider P b acting with respect to 
b-half densities \x~ 1 dxdy\ 1 ^ 2 . After this shift, the relationship between acting on 6-half 
densities and P acting on scattering half-densities is P = xP^x. 

Let VL l J 2 be the bundle of half-densities on M| sc which is spanned by sections of the form 

\f(k,x,y,x',y')^^\^ 2 . 



Let v be a smooth nonvanishing section of this bundle. Since it involves the b-metric g^, 

~ 1/2 

£l b is not the natural bundle near sc. In particular, the kernel of the identity operator on 

M has leading order —n/2 at sc with respect to Q^ 2 |GHlj . We can now define spaces of 
pseudodifferential operators, precisely as in |GHlj : 
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Definition. Let p sc be a boundary defining function for sc. The space ^(M;^ ) is the 
space of half-density kernels K = K\ + K 2 on M| sc satisfying: 

1 ) Psc 2 K\ is supported near Afc jSC , and has an interior conormal singularity of order m at 
Afc :SC; with coefficients whose behavior at the boundary is specified by £; 

2) ps/ 2 -f^2 is polyhomogeneous conormal on M| sc with index family £ , and moreover decays 
to infinite order at bf, lb, and rb. 

The factor of p™/ 2 corrects for the use of b-half densities near sc. Using this definition, we 
can compute as in |GHlj that 

(P + fc 2 )G*f (M;^ /2 ), 

with index sets at sc, 2 at bfo, at zf, 2 at lbo, and 2 at rbo- 
As proven in |GH1] . these spaces satisfy a composition rule: 

Proposition 20. Suppose that A e and B 6 ^ k ' . Then A o B is well defined and 

an element of \&™ +m & ; where 

Gsc = £-sc + J 7 sc, Qzf = (£zf + 3~. z/)0(^ r 6 + J^Wq), Gbf = (£fc/ + ^6/0)0(^0 + -^rbo), 
Glb = (£bf + ^60) 0(^6o + J~zf)i Grbo = (£rb + ^~bf )0(£ z f + J>6o)- 

We therefore expect our parametrix G to be in ^f k for some index family Q. To gain 
more information, we need to start analyzing the model problems. 

4.2. The two-dimensional problem. We begin our analysis of the model problems at the 
face zf. In order to identify the leading-order part of the equation at a boundary hypersurface, 
we need to pick a coordinate to use as a boundary defining function in the interior of that 
face. For all the faces in the lift of {k = 0}, we use k, which is the easiest choice, since it 
commutes with P. Since P + k 2 = xPbX + k 2 , the leading order part of the operator at zf, 
which we call the normal operator, is xPbX, and the model problem is (xPbx)G®j = Id. We 
therefore expect that G Q z j will be (xx')~ l times some right inverse for P&. 

In order to invert Pj,, we use the b-calculus of Melrose [MeJ, identifying zf with the b- 
double space X 2 (x, x') x N y x N y >. The corner zfflbf corresponds to the front face ff in the 
b-double space. An easy calculation, following |GH1] . shows that 

Pb = -{xd x f + (n/2 - l) 2 + A N + W, 

where W is a lower-order term; that is, W vanishes as a b-differential operator at x — 0. 
In fact, Pb is an elliptic b-differential operator, and hence may be inverted by following the 
procedure of Melrose, which is described in |Mej and |Maj . 

The first step in this procedure is to consider the indicial operator, which is the leading 
order part of P& at the front face ff. Using the coordinates (s, x, y, y'), this is 

I ff (P b ) = ~(sd s ) 2 + (n/2 - l) 2 + A N . 

With this terminology, the key theorem is as follows: 

Theorem 21. [MeJ p, is Fredholm as an operator between x s H 2 and x 6 L 2 if and only if 5 
is not an indicial root of Iff(Pb). (See |Mej or |Maj for definitions of x s H 2 and x & L\). 
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In our setting, as in |GH1] . the indicial roots are precisely 

±»i = ± v / (n/2 - l) 2 + \- \ E a(A N ). 

When n > 2, is not an indicial root, and Guillarmou and Hassell show that is not 
only Fredholm but invertible for 5 = 0, and then set G° z j to be (xx')~ l times that inverse. 
However, in our case, n = 2, so z/j = and is an indicial root. So P& is not even 

Fredholm from H 2 to L\. This is precisely why the n = 2 case is not considered in [GHlj. 

4.3. An example: Euclidean space. In order to gain some intuition for the behavior of 
the resolvent near zf in the n = 2 setting, we examine the simplest case, which is M = M. 2 . 
The resolvent on M 2 , acting on scattering half-densities \dz\ = \x~ 3 dxdy\, is 

~H (k\z-zr\) = H (k\y—£\), 

2/T X X 

where H is the Hankel function of order zero. From the asymptotics of the Hankel function, 
we know that Hq(t) decays exponentially as r — > oo, and for small r: 

H (r) ~ -logr + log 2 - 7 + 0(r). 

Using these asymptotics, one can show that the resolvent on R 2 is phg conormal on M 2 SC . 
We are most interested in the leading order behavior near zf. In a neighborhood of zf, we 
have k\z — z'\ < 1 , so this leading order behavior is controlled by the small-r asymptotics 

H (k\z — z'\) ~ — \og(k\z — z'\) + log 2 — 7 = — log k — log \z — z'\ + log 2 — 7. 

Some observations on these asymptotics: 

- As we approach zf, the resolvent increases logarithmically. This is a major difference 
from the n > 3 case studied in [GHlj . in which the resolvent is continuous down to zf. On 
the other hand, the resolvent is continuous down to bf , lb , and rb . 

- The function log \ z — z'\ is the Green's function for the Laplacian on IR 2 . 

These observations suggest that the leading order behavior in the general situation may 
be of the form 

G° zf = -C\ogk + F(x,y,x',y'), 

where F(x, y, x', y') is a right inverse for the operator P with respect to scattering half- 
densities. Moreover, since — Clogfc has logarithmic growth at bfo, lbo, and rbo but the 
resolvent on M 2 does not, we expect that F(x,y,x',y') will have logarithmic growth at those 
faces, with the right coefficient to cancel the logarithmic growth coming from — Clogk. 

4.4. Construction of the initial parametrix. We now construct our parametrix G(k) by 
specifying its leading order behavior at each boundary hypersurface and then checking that 
the models are consistent. 

4.4.1. The diagonal, sc, and bfo. The resolvent has an interior conormal singularity at the 
diagonal {x = x' ,y = y'}. The symbol of P + k 2 is |£| 2 + k 2 , where £ is the dual variable of 
\z — z'\. One can compute that |£| 2 + k 2 is elliptic on M| sc , with leading orders at sc, 2 at 
bfo, and at zf. As in |GH1] . we let the symbol of G(k) be the inverse of |£| 2 + k 2 , in the 
sense of operator composition. This determines the diagonal symbol of G(k) up to symbols 
of order —00, and hence determines G(k) up to operators with smooth Schwarz kernels in 
the interior of M 2 sc . 

At sc, the analysis is identical to that in [GHlj . so we omit some of the details. The key 
point is that sc can be described as a fiber bundle with R n fibers, parametrized by y' G iV 
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and k. The normal operator of P + k 2 is A^n + k 2 , which has a well-defined inverse for k > 0. 
In each fiber, we let 

Gl = (A K 2 + k 2 )~\ 

At bf , we again follow [GHl] exactly. We use the coordinates (k = k/x, k' = k/x',y,y'), 
with k a bdf for bf . Note that these are only good coordinates on the interior of bfo - 
for example, they become degenerate near zf. We then view the interior of bfo as IR + (/t) x 
x N y x N' y . As in |GHlj . the normal operator at bf is 

I b f (kr 2 (P + A; 2 )) = K -\-{Kd K ) 2 + A N + k 2 )k-\ 

Letting Pbf = — (nd K ) 2 + Ajv + k 2 , the model problem is (nPbf K)G^ 2 Q = 5k — n'S—y'. To 

solve it, we separate variables and invert P&/ . For each eigenvalue Xj of Ajv, write Vj = V~\j 
(these are the indicial roots). Let E v . C L 2 (N) be the corresponding eigenspace of A at, and 
let be projection in L 2 (N) onto E u .. Then the inverse of P b f is 

oo 

The only difference between our setting and |GH1] is that we have z/ = as opposed to 
uq > 0. We then set 

G bf = ( KK ')Qbfo- 

We need to check consistency between G° sc and G^; that is, we need to show that they 
agree to leading order in a neighborhood of scflbfo. This proof is the same as in [GHl] ; the 
model problems and formal expressions for G° sc and G^ are identical. We do have u = 0, 
and K (r) has different small-r asymptotics from K v .{r) for Uj > 0; this will be reflected 
in the asymptotics of G^ 2 o near zf. However, since n and k' both approach infinity near 
sc, only the large-r asymptotics are relevant for this consistency check, and K (r) decays 
exponentially just as K v . (r) does. 

Technically, we also need to check consistency between the diagonal symbol and the models 
at bfo and sc. However, this is also the same as in [GHl] ; the models at bfo and sc themselves 
satisfy elliptic pseudodifferential equations given by the leading order part of (P + k 2 )G(k) = 
Id at those faces. As a result, their symbols at the diagonal are determined up to symbols 
of order — oo, and agree up to order — oo with the inverse of |£| 2 + k 2 . 

4.4.2. The leading order term at zf. At zf, the model problem with respect to b-half-densities 
is: 

(xP b x)G° zf = Id. 

Translating our observations in the M = IR 2 case to b-half-densities, we expect 

G% = {xx')-\-C\ogk + F{x,y,x',y')), 

where F(x,y,x',y') is a right inverse for We need to pick the correct right inverse, 
and we expect that it will have logarithmic singularities at all faces. We also need it to be 
consistent with our choice of G^ 2 . To check consistency, we need to show that k°G° z j and 
k~ 2 G b ~f agree to leading order at ff=bfonzf, which is the same as checking if (xx')G° z ^ and 
Qbfo — { KK ')~ 1 C'bf a g ree t° leading order there. 
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First examine the leading order part of Qtf at zf. When x < x', we use the coordinates 
(s, k, x', y, y'), and we have 

oo 

Q bfo = V- 1 I {ks)K {k) +J2^E j Iu j (^)K Uj ( K ). (41) 

i=i 

The boundary defining function for zf is k, so we need to examine the small-K asymptotics. 
For v > we know by standard asymptotics of Bessel functions in |Wa] that 

Io{ r ) ~ 1; K (r) ~ — In r + In 2 — 7; 

U r ) — - — -(-f-KJr) ~ iA!2(l)-". 

V ; r(i/ + 1) 2 ; ' uy 1 2 V 

Here 7 is the Euler-Mascheroni constant. Plugging these asymptotics into (14TI) shows that 
the leading order term in k is: 

7-1 (_ i og K + i og 2 - 7 ) + 

On the other hand, when x > x', we use the coordinates (s' , k' ,x,y,y') and perform the 
same sort of calculations to obtain that the leading order term in k is 



00 



y-i(_log K ' + log2- 7 ) + ]rn^ 



3=1 1 



The x < x 1 and x > x' cases may be combined; we see that the leading order term of Qbf 
at zf is 

V- 1 (- log k + log 2 - 7 + log x' + x s <i log s) + J2 n E Vj ' ■ (42) 



3=1 J 



We need to construct G° z j- so that (xx')G° z j- has this leading order at bfo- 

To construct G° z ^ we must first find the correct right inverse for p,. Fix S with < 5 < V\\ 
then P& is Fredholm from x~ s H 2 to x~ s L 2 by Theorem [2TJ Following the usual b-calculus 
construction in [MeJ and |Maj . we obtain a generalized inverse Qt S . We claim: 

Lemma 22. P5 zs surjective onto x~ s L 2 . 

The lemma implies that Q^ 6 is an exact right inverse for p,. 

Proof. By taking adjoints, the lemma is equivalent to the statement that P& is injective on 
x 5 L 2 . Suppose that u\dgb\ 1 ^ 2 is in x 5 L 2 and satisfies P b u = 0. By regularity of solutions to 
b-elliptic equations, u is phg conormal on M near x = 0; since u G x s Lf, it decays to at 
least order 6 at 2 = 0. On the other hand, since P = xPji and Idgt] 1 ^ 2 = x\dg\ l l 2 , u\dg\ l l 2 
is in the kernel of P = Am, and hence Aj^u = 0. By the maximum principle, u — 0, which 
completes the proof of the lemma. □ 

The correct right inverse will be a slight modification of Q^ S ■ In order to check consistency, 
we need to understand the structure of near the front face ff=bfonzf. This structure is 
described in detail in |Me] and |Maj . In particular, the leading order of Q^ s at ff is precisely 
the indicial operator Iff(Q^ S ), which satisfies the equation 

- (sd s ) 2 + A N )I ff (Q b ) = 6(s = l,y = y'). (43) 
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Moreover, from [Mej and [Ma], Iff(Qb~ S ) nas polyhomogeneous expansions at s = and 
s = oo, with leading order terms at worst s~ 5 at each end; that is, a small amount of growth 
is allowed at s — 0, and a small amount of decay is required at s = oo. 

We now separate variables and solve (1431) directly. For each j > 1, s ±Uj span the kernel 
of — (sds) 2 + Uj. Therefore, the solutions corresponding to E v . are combinations of s v i and 
s~ U] away from s — 1. By the requirements at s = and s = oo, our solution is a multiple 
of s~ Uj for s > 1 and of s Vj for s < 1. Using the matching conditions at s = 1 arising from 
the delta function singularity, the solution on the eigenspace E v . is 



n 



e 



^ | logs | 



E v . 



2uj 

We have to consider z/ = separately; the kernel of — (sd s ) 2 is spanned by 1 and logs. 
Because we require decay at s = oo, the solution for s > 1 must be zero. Then the matching 
conditions at s — 1 imply that the solution is logs for s < 1. Since projection onto E is 
simply V -1 , where V = Vol(N), the zero-eigenspace solution is Vx{ s <i} logs. Therefore the 
leading order part of Q7 S at ff is 



Vj | logs | 



iffiQt 5 ) = v^x^niogs + ^Tn^^-— . (44) 

3=1 J 

Now compare (j44p with (1421) . We see immediately that if we let 

= (xxO^^-U" 1 log A; + Q~ s + V~ x logx' + V 1 (log2 - 7)), 

then and G^f are consistent. In addition, since any function independent of x and y 
is in the kernel of Pb, G° z j- solves the model problem (xPbx)G Q z j = Id at zf. Moreover, the 
diagonal symbol is consistent with for the same reason that it is consistent with G° sc and 

G bf - 

4.4.3. The model terms at rb^. Finally, we need to specify the leading-order behavior of the 
parametrix at rbo; in fact, we need to specify some lower-order terms as well. We use the 
coordinates (x, y, k', y'); the k' = face is rboflzf and the x = face is rboflbfo. There will be 
a term G^ bo 1 for each Uj in [0, 1). The model problem near this face, with A; as a boundary 

defining function, is (xPbx)u = 0, so we need Pb(xG rb U3 ) = for each Uj 6 [0, 1). 
First we focus on the model of order —1. We let 

and claim that this is consistent with G° z ^ and G b f Q . 

To check consistency with G z j, we need to show that the leading order of G° z f agrees with 
the leading order of k~ 1 G~ b 1 Q at zfflrbo. Recall that at rbo, which corresponds to s = 00, 
Q b s decays to a positive order. So (xx'^Q^ 5 has leading order greater than —1 at rb ; 
therefore, the leading order part of G° z j at rbo is precisely 

{xx')- l {V- l {- log*/ + log 2 - 7)). 

But by Bessel function asymptotics, 

k ~ lG rb\ = V-\xx')~ 1 K (k') ~ (xx')-V- 1 (-logfi: , + log2-7). 

Therefore G° z j and G~ bo are consistent. 
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We must also check consistency of G rbo with G h f. Near rb 



i=i 

We are only interested in the order —1 part of this term. Since x' and k both vanish to first 
order at rbo, all the j > terms have leading order —1 + Uj at rbo- Since Io(0) = 1, the 
order —1 part of k~ 2 G b 2 Q at rbo is precisely 

{xx')- 1 V^K (k') = k^iV^x^K'K^x')) = k^G^l- 

We conclude that G~ bo is consistent with the models at zf and bf . 

We also need to specify some lower order terms at bfo; for this we precisely follow Section 
4 of |GH1]. At zf, they need to match with the asymptotics of Q b s , and at bfo, they need to 
match with the higher order Bessel functions. Both of these involve only the nonzero indicial 
roots, so the terms and arguments are identical to [GHlj . In particular, for any < Uj < 1 
in the indicial set, we let 

Q u s -l _ i iJK Vj (K?) 

^rb - x r(v j )2"s-i Vi{ ' y} ' 

where v Uj (z, y') is in the kernel of P& with asymptotic 

v Vj {x,y,y') = (2u j )- 1 U Ej x-^ +0(x^- 1 \ogx). 

The function v Uj is there to match with the asymptotics of Q b at rbo, as i n Section 4 of 
[GHlj . In fact, these models are consistent with our models at bf and at zf by precisely the 
same argument as in [GHlj; we will not repeat it here. 

4.5. The final parametrix and resolvent. We have now constructed models at sc, bfo, zf, 
and rbo which are consistent with each other and also with the diagonal symbol. Moreover, 
all the models decay to infinite order as we approach lb, rb, or bf. Therefore, we specify 
our parametrix G{k) to be any pseudodifferential operator in with kernel having the 

specified diagonal symbol and specified leading-order terms at sc, bf , zf, and rb . The 
consistency we checked guarantees that such an operator exists. The behavior of the kernel 
of G(k) at lbo may be freely chosen as long as the leading-order term is order —1 and it 
matches with our models at zf and bfo; a term of order —1 will, however, be required. 

Now let E(k) = (P + k 2 )G(k) — Id. Since G(k) has diagonal symbol equal to the inverse 
of the symbol of P + k 2 , the Schwarz kernel of E(k) is smooth on the interior of sc . 
Moreover, since P + k 2 is a differential operator, the Schwarz kernel of E(k) is phg conormal 
on M 2 . 

w,sc 

- At lb, rb, and bf, the Schwarz kernel of G(k) vanishes to infinite order along with all 
derivatives, so the same is true of E(k). 

- At sc, G° sc solves the model problem, so E(k) has positive leading order at sc. 

- At bfo, G{k) has order —2, but G^f solves the model problem, and moreover P + k 2 
vanishes to second order. Therefore E(k) has positive leading order at bfo- 

- At zf, G° z j solves the model problem, so E{k) has positive leading order. 

- At lbo, G(k) has order —1. The variables k and x both vanish at lbo, so k 2 G(k) has order 
1 and xG(k) has order 0. Since P b is a b-differential operator, Pf,(xG(k)) also has order 0, 
and hence (P + k 2 )G(k) = {xP^x + k 2 )G(k) has order 1. Since Id is supported away from 
lb , E(k) decays to at least order 1 at lb . 
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- At rbo, G(k) has order —1, but all the terms G Vj ~ 1 for Uj G [0, 1) solve the model problem. 
Therefore, the error E{k) has leading order at worst 0. 

To summarize, if £ is the index set for E(k), we have shown: 

£ sc > 0, £ z f > 0, £bf > 0, £ib > 1, £ r b > 0, £ib = £ r b = £bf = 0. 

Now we iterate away the error. By Proposition [201 E(k) 2 vanishes to positive order at all 
faces of M| sc ; suppose that the order of vanishing at each face is greater than e > 0. Again 
applying Proposition [20| we see that for each JVgN, the order of vanishing of E(k) 2N and 
E(k) 2N+1 at each face of M| sc is greater than N5. Therefore the Neumann series 

oo 

(id+Eik))- 1 = Y, E ( k ) 1 

may be summed asymptotically, and the sum defines an element of for some index 

family £ . 

Finally, let R(k) = G(k)(Id + E(k))~ l ] we see that (P + k 2 )R(k) = Id. Since P + k 2 is 
invertible for all positive k, its only right inverse is the resolvent. We conclude that R(k) is 
in fact the resolvent, and it is an element of %°°' n (M, S\ 1/2 ) for some index family 1Z. 

In order to prove Theorem El we need to perform this construction for any angle 9 G 
(— 7r, 7r), not just for 9 = 0. However, as Guillarmou and Hassell claim in iCillTI. the con- 
struction is essentially unchanged. Indeed, we just use 

e i9/2 k 

as our boundary defining 

function for the k = faces instead of k, and correspondingly change the model at sc from 
(A R 2 + k 2 )^ 1 to (A R 2 + e ld k 2 )~ l . The construction is then precisely analogous to the 9 = 
case; by construction, the index sets are independent of 9. Moreover, by the continuity of 
the resolvent outside the spectrum (also by construction), all the dependence on 9 is smooth. 
This completes the proof. 

4.6. Leading orders of the resolvent. Since R(k) = G(k) -G(k)E(k) + G(k)E(k) 2 - . . ., 
we can obtain some information about the leading orders of R(k) at each face. For n > 3, it 
is shown in [GHlj that the leading orders of R(k) are the same as those of G(k); we claim 
that the same is true when n = 2. 

When n = 2, G{k) has leading orders —1 at lbo and rbo, order at sc, order -2 at bfo, and 
logarithmic growth at zf. E(k) has non-negative leading orders at all faces, and it is easy 
to use Proposition [20] to show that the leading orders of G{k)E{k) are no worse than those 
of G(k). Similarly, it may be shown that the leading orders of G(k)E(k) 1 are no worse than 
those of G(k). Since G(k) is fixed and Id — E(k) + E(k) 2 — ... is asymptotically summable, 
the series G(k) — G(k)E(k) + . . . is also asymptotically summable; therefore, the leading 
orders of R(k) are no worse than those of G(k). The leading order terms themselves may be 
affected by G(k)E(k), but the orders are not. 

To summarize, when n = 2, the leading orders of the exact resolvent R{k) are no worse 
than at sc, -2 at bfo, logarithmic at zf, and —1 at lbo an d rbo- When n > 3, the orders are 
at worst at sc, -2 at bf , at zf, and | — 2 at lbo and rbo, as in [GHlj. However, these orders 
are for the resolvent acting on b-half-densities, rather than the more natural scattering half- 
densities. Switching to scattering half-densities requires an order shift, adding n/2 at each 
of {x = 0} and {x' = 0}. So we need to add n/2 to the orders at lb and rb , and n at bf . 

The order at sc remains unchanged, because the extra factor of pl^ in the definition of the 
calculus already incorporates the shift. So: viewing the resolvent as a scattering half-density 
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\dgdg'\ l l 2 acting on scattering half-densities for each k, or equivalently as a function acting 
on functions on M by integration against dg, it has leading orders given by: 

• at sc and n — 2 at bf , rbo, and lbo: 

• r z f at zf, where r z f = if n > 3 and r z f = (0, 1) (that is, leading order behavior of 
log p zf ) if n = 2. 

Appendix A. Construction of the short-time heat kernel 

Albin has created a framework for the construction of the heat kernel on an asymptotically 
conic manifold; essentially all of the hard work involved in this construction has already 
been done in [Alb] . To complete the construction and prove Theorem [31 all we need to do 
is create an initial parametrix for the heat kernel. This construction is the content of this 
short appendix and is based on Section 5 of [Alb] . in which Albin constructs the heat kernel 
on an edge manifold. 

The space in Theorem [31 which we call S, is obtained by taking the manifold M 2 C x [0, T) t 
and then blowing up the t = diagonal. We call the scattering face of M 2 C x [0, T) t sf and 
call the front face at t = ff . The heat operator is precisely d t + Am- Our goal is to create 
a parametrix which, to first order, solves the normal equations at sf and ff. 

First analyze the situation at sf . As first discussed in |Me3j and elaborated upon in [GHlj , 
sf has a Euclidean structure, parametrized by y £ iV and t £ [0, T). By the same analysis as 
in |GHlj , the normal operator at sf is precisely d t + A K n . We then simply let the model at 
sf be the Euclidean heat kernel, H Rn . This is analogous to the construction in Section 5 of 
[Albj for the edge setting, in which the model is the heat kernel on hyperbolic space times 
the heat kernel on the fiber. 

The analysis at ff is also standard, since ff corresponds to the short-time regime on the in- 
terior of the manifold, where the heat kernel asymptotics are local. We know that d(z, z') / yt 
is a good coordinate along ff, zero at the spatial diagonal and increasing to infinity as we 
approach the original t = face. We therefore let the leading-order model at ff be 

(rf(*,z')) 2 

The choice of model at ff is again based on the Euclidean heat kernel, and is precisely the 
same as the choice of model in the edge setting |Alb] . 

Each model vanishes to infinite order as we approach all boundary hypersurfaces other than 
sf and ff. Moreover, the models are consistent, as the leading orders of each are precisely the 
Euclidean heat kernel at sf PI ff. We may therefore pick a pseudodifferential operator whose 
Schwarz kernel agrees with our models to leading order at sf and ff, and decays to infinite 
order at all other boundary faces. In Section 4 of [Alb] . Albin proves a composition rule for 
time-dependent pseudodifferential operators whose kernels are polyhomogeneous conormal 
on S - our setting is the 'scattering' setting, which is included in his analysis. We then 
use this composition rule and an iteration argument, precisely as in Section 5 of [Alb] . to 
construct the heat kernel as a polyhomogeneous conormal distribution on 5*. This completes 
the proof of Theorem [31 
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